N-order bright and dark rogue waves in a resonant erbium-doped fibre system 
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The rogue waves in a resonant erbium-doped fibre system governed by a coupled system of the 
nonlinear Schrodinger equation and the Maxwell-Bloch equation (NLS-MB equations) are given 
explicitly by a Taylor series expansion about the breather solutions of the normalized slowly varying 
amplitude of the complex field envelope E, polarization p and population inversion r;. The n- 
order breather solutions of the three fields are constructed using Darboux transformation (DT) by 
assuming periodic seed solutions. What is more, the n-order rogue waves are given by determinant 
forms with n + 3 free parameters. Furthermore, the possible connection between our rouge waves 
and the generation of supercontinuum generation is discussed. 
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I. INTRODUCTION 

In recent years, long haul optical communication 
through fibers has attracted considerable interest in re- 
search activities among scientists all over the world. Es- 
pecially, it has been demonstrated that the soliton-type 
pulse propagation will play a vital role in the ultra fast 
communication systems. They are considered to be the 
futuristic tools in achieving low-loss, cost-effective high 
speed communication throughout the world. Soliton- 
type pulse propagation through nonlinear optical fibers 
is realized by means of the exact counterbalance between 
the major constraints of the fiber, viz., group velocity dis- 
persion (linear effect) which broadens the pulse and the 
self-phase modulation (nonlinear effect) which contracts 
the pulse. The propagation of optical pulses through a 
nonlinear fiber in the picosecond regime is described by 
the well-known nonlinear Schrodinger (NLS) equation, 
which was first proposed by Hasegawa and Tappert in 
19731]. 

To make the soliton based communication systems 
highly competitive, reliable and economical when com- 
pared to the conventional systems, attenuation in a fiber 
must be compensated. A different type of optical soliton 
is associated with the self-induced transparency (SIT) 
effect in resonant absorbers. The soliton pulse prop- 
agation in an erbium doped fiber amplifier utilizes the 
SIT phenomena, first discovered by McCall and Hahn 2] . 
In 1967, McCall and Hahn proposed a new type of op- 
tical soliton in a two level resonant system. Above a 
well-defined threshold intensity, short resonant pulses of 
a given duration will propagate through a normally ab- 
sorbing medium with anomalously low attenuation. This 
happens when the pulse width is short, compared to the 
relaxation times in the medium and the pulse centre fre- 
quency is in resonance with a two-level absorbing transi- 
tion. After a few classical absorption lengths, the pulse 
achieves a steady state in which its width, energy and 
shape remain constant. The pulse velocity has greatly 
reduced from the normal velocity of light in such media. 
With these properties, the pulse propagation of this type 



is named " self- induced transparency" (SIT) soliton and 
frequently described by the Maxwell-Bloch (MB) equa- 
tions. They are 

Pt =i uop - fqv, 

Vt^fiqp* + q*p), (1) 

Here, E and p are complex variables, 77 is a real variable 
and Wo is a real constant and / is the character describ- 
ing the interaction between the resonant atoms and the 
optical field. The symbol * denotes the complex conju- 
gate. These equations can be extended to the case of 
fiber amplifiers. When Er is doped with the core of the 
optical fibres, then the nonlinear wave propagation can 
have both the effects due to silica and Er impurities. Er 
impurities give SIT effect to the optical pulse, whereas 
the silica material gives the NLS soliton effect. So if we 
consider these effects for a large width pulse, then the 
system dynamics will be governed by the coupled system 
of the NLS equation and the MB equation (NLS-MB sys- 
tem). Considering the Erbium doped in nonlinear silica 
waveguides, for the first time, the combined NLS-MB 
system was proposed by Maimistov and Manykin [3, Hjin 
1983. They have also constructed the Lax pair and used 
the inverse scattering transform technique for the gener- 
ation of soliton solution. The NLS-MB equations read 
as[^-^5,] 

Et^i[^E.,, + \E\''E]+2p, 
Px =2i loqp + 2Eri, 

in.,^-{Ep* +E*p), (2) 

The above equation have also been reduced through the 
Painleve analysis Q. Further, Kakei and Satsuma Q 
also reported the Lax pair and the multi-soliton solu- 
tion of the NLS-MB equations. The integrability aspects 
of NLS-MB system with variable dispersion, the study 
of propagation of optical solitons in coupled NLS-MB, 
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and random nonuniform-doped media have been reported 
earlier wherein the spectral parameter was kept constant. 
The coexistence of NLS soliton and SIT soliton has al- 
ready been confirmed experimentally [lOi [flj . The propa- 
gation and switching of SIT in nonlinear directional cou- 
plers with two-level atom nonlinearity has been recently 
investigated numerically by retaining the transverse de- 
pendence of the optical field and atomic variable. Recent 
experiments by Nakazawa et al., have confirmed guided 
wave SIT soliton formation and propagation by employ- 
ing a few meters of erbium doped fibre [8l-[T]| . 

Recently, considering all higher order effects in the 
propagation of femtosecond pulses, the coupled Hirota 
and Maxwell-Bloch (CH-MB) equations have been pro- 
posed and analyzed for soliton solutions . Some gen- 
eralization of NLS-MB equations, for instance, the CH- 
MB equations and the NLS-MB equations with variable 
dispersion and nonlinear effects are discussed Ta'-'Tl. 
The single soliton and the single breather solutions 16j of 
the NLS-MB equations are given by the Darboux trans- 
formation (DT)[l3, [13 • Soliton solution for the gener- 
alized coupled variable coefficient NLS-MB system was 
also investigated by the DT [I9| and the Hirota method 

In recent years, in addition to solitons in different op- 
tical systems, the study of rogue waves have also at- 
tracted considerable interest because of their potential 
applications in different branches of physics including 
oceanographv|2l|-|23|. which occurs due to either mod- 
ulation instability |25l431| . or random initial condition 
[2^ [3^ . The first order rogue wave is most likely to ap- 
pear as a single peak hump with two caves in a plane 
with a nonzero boundary. One of the possible gen- 
erating mechanisms for rogue waves is the creation of 
breathers which can be realized by modulation instabil- 
ity. Then, larger rogue waves can build up when two or 
more breathers collide themselves [33l - l39j . Recently, more 
general higher-order rogue waves were obtained such as 
showing that these general N-th order rogue waves con- 
tain (A^ — 1) free irreducible complex parameters 
Rogue waves can also be observed in space plasmas 
fisi and optics when propagating high power optical radi- 
ation through photonic crystal fibers [46l-|48|. Consider- 
ing all higher order effects in the propagation of femtosec- 
ond pulses, rogue waves can also be observed in a system 
modeled by the Hirota equation [49l - l5l| . Furthermore, 
rogue waves have not only been observed in continuous 
media but have also been reported in discrete systems, 
such as the systems of the well-known Ablowitz-Ladik 
(A-L) equation (52j . 

Though the rogue waves have been reported in differ- 
ent branches of physics where the system dynamics is 
governed by single equation, to the best our knowledge, 
they have been observed and reported very little in 
the coupled systems. For example. Rogue waves of the 
coupled NLS were constructed in the literatures 's^-'FE'] . 
Very recently, new kinds of matter rogue waves (56] have 
been reported in F = 1 spinor Bose-Einstein condensate 



37 . 

4C- 



system controlled by three component NLS equation. 
In experiment, the rogue waves in a multistable system 
is revealed by experiments with an erbium-doped 
fiber laser driven by harmonic pump modulation. 
So, it is our prime interest to analyze the possibility of 
rogue waves in coupled systems, such as NLS-MB system. 

It is well known that the dark soliton [s^ of the 
defocusing nonlinear Schrodinger (NLS) equation is 
essentially different from that of bright soliton. For 
the past two decades or so, intensive research have 
been carried by several groups about theoretical and 
experimental aspects of dark bright solitons, it is quite 
natural to ask a question: is there any possibility of 
observing dark rouge wave in soliton equations? In 
general, the first order dark rogue wave has one down 
dominant peak and two small lumps. Because of the 
singularity H!] of the solution for the de-focusing 
NLS equation generated by using DT, we cannot 
get dark rogue wave of the de-focusing by this way. 
Fortunately, we have obtained dark and bright rogue 
waves 59J of the NLS-MB equations from a Taylor series 
expansion of the first order breather solutions, which 
are generated from a periodic seed by the DT. But we 
did not provide a detailed analysis on their dynamical 
evolution and higher order rogue waves. The aim of this 
paper is twofold. Firstly, the determinant representation 
of the n-fold DT of the NLS-MB equations are similar 
to the case of NLS equation DT[60]. Secondly, the 
rogue waves of the three optical fields are constructed 
by determinant forms. It should be noted that the rogue 
waves of the fields p and rj are dark. Furthermore, the 
connection between our rogue waves and the generation 
of supercontimmm generation will be discussed. 

The organization of this paper is as follows. In section 
2, the determinant representation of the n-fold DT and 
formulae of iS'"! , p^"! and t^I"' are expressed by eigenfunc- 
tions of spectral problem. In section 3, a Taylor series 
expansion about the breather solutions are generated by 
n-fold DT from a periodic seed solution with a constant 
amplitude to construct the bright and dark rogue waves. 
What is more, the n-order rogue waves are given by de- 
terminant forms with n + 3 free parameters. Finally, we 
conclude the results in section 4. 



II. DARBOUX TRANSFORMATION 

The linear spectraLproblem of the NLS-MB equations 
can be expressed as 



(3) 
(4) 



where 
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and A is the complex eigenvalue parameter. 
It is easy to prove that the spectral problem ^ and @ 
are transformed to 



The elements of one-fold DT [T^ are parameterized by 
the eigenfunction /fc associated with A^ as 
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under a gauge transformation 

v|/[i] = r vj/. 



(7) 



Here, T is a 2 x 2 matrix, which is determined by the 
cross differentiating (O and ([6]), 

f/w* -v^w. + T/W] = r([/t- 14 + [f/,V])T-i. (8) 

This implies that, in order to make eq.® and eq.(|3]) 
invariant under the transformation it is crucial to 
search a matrix T such that J/'^l and y'^' have the same 
forms as U and V. At the same time the old potential 
(or seed solutions) {E, p, rj) in spectral matrixes U and 
V are mapped into new potentials (or new solutions) 
ii-y terms of new spectral matrixes [/[^l 

and yW. 

2.1 One-fold Darboux transformation of NLS-MB 
equations 

In order to be self-contained, we shall recall the 
one-fold DTfTe'] of NLS-MB equations. Considering the 
application of the representation for the n-fold DT by 
means of the determinant [60j of eigenfunctions with 
different eigenvalues in the following context, we need to 

/fei 

fk2 

associated with an eigenvalue Afc,and Xk = A™ if fc = m, 
where k — 1,2,3, ...,2n but Afc ^ X. Additionally, 
the eigenfunctions for distinct eigenvalues are linearly 
independent, (i.e.) fk and fm are linearly independent 
if fc 7^ m. 



introduce 2n eigenfunctions by fk = fk{Xk) = 
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With the transformed potentials. 
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The resulting new solutions of £'[^1 , p^^l and rp-"^ are given 

by 
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.[1] = 



1 



det{Ti) 



(-277(ri)ii(Ti)i2 +p*(ri)i2(ri)i2 -p(ri)ii(Ti)ii)|A= 



det{Ti 



-(?7((Ti)n(Ti)22 + (Ti)i2(Ti)2l) -p*(ri)i2(Ti)22 +p(T^l)ll(Tl)2l)U=^, 



(11) 

(12) 
(13) 



and the new eigenfunction fj^^ of Xk corresponding to the 
new potentials is 
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In order to satisfy the constraints of 5' and VLi in [T6| . 
set 
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2.2 n-fold Darboux transformation for NLS-MB equa- 
tions 

In this subsection, our prime aim is to estabhsh the 
determinant representation of the n-fold DT for NLS- 
MB equations as we have done for the case of NLS 
eauation(60j. According to the form of Ti in eq.®, 
the n-fold DT should be of the form T„ = Tn{X) = 

A"-/ + hX"-^ + t2X"^^ H h t„-iA + tn, where U are 

2x2 matrices, i = l,2,---n. T„ which leads to the 
determinant representation of r„ by means of its kernel. 
Specifically, from algebraic equations, 



coefficients ti, I = 1,2, ... ,n are solved by Cramer's rule. 
Thus we obtain the determinant representation of the T„ . 
Theorem 1. The n-fold DT of the NLS-MB equations 
is r„ = r„(A) = A"/ + tiA"-i + t2A"-2 + • • • + tn-iX + 
i„,where U are 2x2 matrices, i — 1,2, ■ ■ ■ n. The final 
form of Tn{X) has the form. 



\W2n\ \W2n\ 

Tn ~ Tn{X; Ai, A2, • • • , A2n-1, A2n) = 

{Tn)21 {Tn)22 

\ \W2n\ \W2n\ J 

(16) 

Here / is a unit matrix, and expression for ti and com- 
ponents of Tn are given in Appendix I 

It is easy to construct a simple form of the determinant 
ofT„ 



det{Tn) = (A - Ai)(A - A2) • • • (A - A2„-i)(A - A2„) 



Next, we consider the transformed new solutions 
(i?["'l,p["'l, T?!"!) of NLS-MB equations corresponding to 
the n-fold DT. 

Corollary] l.For the n-fold DT, the transformed poten- 
tials are 
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n 
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which leads to the new solutions i?^"' and r;!"! of the 
form 



5 
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Note that 

A2fc = -A2fe_i, /2fc = ( Z^'""^^) (21) 

V J2A:-11 / 

in order to satisfy the constraints of DT. 



ated by assuming a periodic seed solution. Then we can 
construct the exphcit bright and dark rogue waves of the 
NLS-MB equations through a Taylor series expansion of 
the breather solutions. 



III. N-ORDER BRIGHT AND DARK ROGUE 
WAVES GENERATED BY N-ORDER 
BREATHER SOLUTIONS 

By using the results of DT discussed above, breather 
solutions of p and 77 of NLS-MB equations are gener- 



Substituting, E — (iexp[i/9],p — ifE,rj — 1 into the 
spectral problem eq.(|31) and eq.(|ll), and using the method 
of separation of variables and the superposition principle, 
the eigenfunction /2fc-i associated with A2fc_i is given by 



(hk-ll{x,t,\2k-l\JCiw{x,t,\2k-l)[l,2k - I] - C2W*{x,t,-\l^_^)[2,2k - l\ 
l/2fc-i2(a:, t, \2k-1)) lCitz7(x, t, A2fc-i)[2, 2fc - 1] + C2VJ*{x, t, -A*^,_i)[l, 2k - 1]) ' 

Here 

/m{x,t,\2k-i)[l,2k- l]\f dexp[^p + ic{X2k-i)] \ 

\^ix,t,\2k-i)[2,2k-l])-\[t (ci(A2fc_i) + I) -"A2fc-i] cxp[-|p + zc(A2fc-i)y ' 

, , , / ti7(a;,t,A2fc-i)[l,2fc- 1] \ 

U(^'*'^2fc-i)[2,2fc-i] 

I 



Note that -07(2;, A2fc-i) is the basic solution of the spec- tral problem eq.® and eq.(|l]). Here a,b,d,t,z e M, 



6 



tions 



c{\2k- 


-l) 


=Cl(A2fe- 


-i)x 


Cl(A2fc- 


-l) 


V 


.ib 

^ 9 

z 


C2(A2fc- 


-l) 




-1 ~ 




P 


=at + b 


X, 




f 


1, 


2 " 



A2fc-1 — »Wo 



Cl(A2fe-l), 



and the {—b + 2ujo)f = 2 is used for /. 

3.1 The n-order breather solutfons of NLS-MB equa- 



For simplicity, under the condition Ci = C2 = 1, let 
A2fe-i a2k-i + ^^'^^ that Im(^ - A2fe-i) = and 



ci(A2fc-i) = '\J(P - Kc^{\2k-i) e R,and using eq. 
then substituting eigenfunctions eq. into eq. (|16l) . 

we obtain the determinant representation of DT in the 
form which is discussed in Appendix II. 



Using eq. (|A.ip into eq. (|18ll9l20p with the choice in 
eq. ()A.2p . we can construct and Tyl"'. For brevity, 

in the following, we are giving only an explicit expression 

of with specific parameters d= \,b — 2,uJo — -. 



J 



=E + 4ai— exp(i(-5i + 2a;)), 



vz = —ai cos(2wi) + (2 cos(t«i) + 2\J 1 — a\ sin(?i;i)Q;i + 2ai^ cos(it;i)) cosh(w2) 



+ (2iai2 sin(wi) — 2isin(u)i) — 2iQ;i cos(wi)-\/l — ai^) sinh(w2) — 2q;i — ~ ai^ sin(2wi) 
+ — sinh(2w2) — ai cosh(2u'2), 

W4 = 2qi2cos(2wi) — 2(4ai cos(w;i) + 2sin(wi)-\/l — ai^) cosh(w2) 
— 2(— 1 — ai^ — cosh(2w2) — vl — Oi^ai sin(2wi)). 
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wi = —2\J\ — ai^x + — -v/l - ai^t, W2 = —\[\^^oi{^t. 

5 5 



The dynamical evolution of jp'^lp and ryl^l for 

the parametric choice c? = 1,6 = 2,wo — 2'*^! ~ ^'"^ 
are plotted in the Figures 1-3, which confirms the direct 
verification of the periodic as well as decaying properties 
of typical breather solutions. The breather of E^^^ is 
almost same as that of the NLS equation, which has one 
upper peak and two caves in each periodic unit. On the 
other hand, we observe that there are two new kinds of 
breathers for pl^l and r^^^. It is interesting to note that 
new breather pl^l admit one upper ring and three down 
peaks in each periodic unit. Whereas the new breather 
ryl^l has two lumps and one down peak in each periodic 
unit. Moreover, these two new breathers can be called as 
dark breathers because the down amplitude is dominant 
in both the cases. The above discussed new properties 
are clearly seen in Figures 1-3. 



3.2 The first-order rogue waves generated by first-order 
breather solutions above 

Similarly, under the condition Ci = C2 = 1, sub- 
stituting eigenfunctions ea. (l22t into eqs. (ll8ll9l20l) with 

Ai = ai + i-, by assuming a\ — > d[d > 0), E^^\p^^^and 



r?t^] become rational solutions {M^\ p^^\ fj^^^ in the form 

t 



of rogue waves [59|. When x ^ 00 



above expressions, after some manupulations, we find 
^ d^, |pW|2 ^ and ryW ^ 1. In 



00 in the 



addition to the above conditions, from I-E'^'I^ = and 
= 0, we also observe that the maximum ampli- 
tude of jijt^lp occurs at t = and x — and is equal 
to 9d^, and the minimum amplitude of j^I^lp occurs at 
a/3 

t = and X — and is equal to 0. By using sim- 

ilar procedure discussed above, we can also obtain the 
extreme value of jp'^lp and t^'^I. 

Figure 4 is plotted for the rogue wave ^ with spe- 
cific parameters d = 1, 6 = 2, = — . From figure 4, we 

infer the following interesting results: 1) the ji^I^lp 1 
by assuming x ^ 00, i — > 00 which gives the asymptotic 
plane; 2) The maximum amplitude of jiJ^^lp occurs at 
t — and X — and is equal to 9, and the minimum 

amplitude of \E^^^ p occurs at t = and x = and is 

equal to 0. As the general expression of the extreme val- 
ues of 1 2 and fj^^^ are quite complicated in nature, for 
simplicity, we only discuss these solutions under certain 
choice of parameters. 

Figure 5 is plotted for the rogue wave p on {x — t) 
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plane with the above parameters. Like in the earher 
case, here also we observe the following salient features: 
1) the height of the asymptotical plane is 4 because 
_^ 4^ when x — >■ oo, t — >■ oo; 2) The maximum 
amplitude of occurs in the form of ring curve on 

{x — t) plane defined by 
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occurs at tour points { (i = — , x = 
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equal to 0; 3) the extreme value of the amplitude 

4 

occurs at t = and x — and is equal to — . We also 

Zo 

observe that the middle down peak in Figure 5 has two 
sub-peaks. Due to the direction of the observation of 
the figure, these two close sub-peaks are not clearly dis- 
tinguished from the figure, we just find three down peaks. 

Figure 6 is plotted for the rogue wave 77^^! with 
specific parameters as in fig. 4. From the Figure, we 
observe the following new results: 1) the height of the 
asymptotical plane is 1 because t/'^I ^ 1 by letting 
a; — > 00, f — > cx); 2) the maximum amplitude of ryl^l 

occurs at two pomts {{t — H — — ,x = H — j 

X A su 5 + %/5 19 + 9^ ^ 

I and i (i = — — ,x = — ) I and is equal 

52 52 

to \/5, and the minimum amplitude of fj^^^ occurs at 

• , su +19-9^5 
two points {[t = ,x — ) I and 



{it = 



V5 



52 



52 

-19 + 9V^ 
52 ' 



52 



} and is equal to —Vb; 



3) the extreme value of the amplitude rj^^^ occurs at i = 

and a; = and is equal to . Like in Figure. 5, here 

5 

also we observe that the down peak in Figure 6 has two 
sub-peaks. 



3.3 The higher order rogue waves and their determi- 
nant forms 

In order to emphasize the richness of the higher order 
rogue waves, we can modify Ci and C2 in the equation 
as the following: 



fe-i , 
)} and is C,=Ko + exp{zci{X2k-i)J2jjiX2k-i-id + i-)y) 



fc-i 



C2 = Ko + exp(-ici(A2fe-i) ^ J,-(A2fe-i - id + i-)y) 

3=0 

(23) 

h 

Here Kq, Jj e C. Note that A2/C-1 = d 
point of Ci(A2fe-i). 



z— is the zero 
2 



Based on the section 3.2, higher order rogue waves can 
be constructed by the breather solutions. In other words, 

let A2fc-i — > d + i— in n-order breather solutions, n-order 

rogue waves can be given. Generally, in comparison to 
the method of limiting the breather solutions, the method 
of making rational eigenfunction below may be more di- 
rect and the rogue wave can be shown by determinant 
forms. 

Substituting eg. ([25)1 into eas. ip^ . by assuming 

A2fe-i d + i — , eigenfunction /2/c-i associated with 
A2fc-i become rational eigenfunction fr as follows. 



ii2dKo + 2d)x+2{- 



2i 



-+i 



d - b){Ko + l)dt + 2Jod + Ka + l)\/2dexp(i^) \ 



(& — 2ujQ)[d + —lb — iujo) 
-{2dKo + 2d)x-2{ '^-^ +id - b){Ko + l)dt - 2Jod + Kq + l)V2dexp{-K] 



(b — 2uJo)id + —ib — iujq) 



J 



1 , 8 + b^ - 2b'^ujo - 2d'^b + M'^ujo , , 
^ - 2*( ' + 



(24) 



Substituting eigenfunctions ea.([M| into eas. (|llll2ll3p . 

we can get the first order rogue waves {E^^\p^'^\fp-^} 
in the form of determinant. The dynamical evolu- 
tion of |£'[^lp, Ip'^'P and t/'^' for the parametric choice 



d, 6, Wo, Kq, Jq are respectively similar to the Figures 4-6, 
but we can control the position of the first-order rogue 
waves by choosing the parameters Kq and Jq. 

Theorem 2. For the n-fold DT, the n-order rogue waves 
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and of the form 



P'"' = 1 j-trf. N (-2?7(7V„)ii(7V„)i2 +P*('7V„)l2(7Vn)l2 -p('f'™)ll(7Vn)ii)|A=ic^o, 

det{l rn) 



(25) 
(26) 

(27) 



r 



The final form of r„i(A) is given in Appendix III. 
Case 1). When n — 2, substituting eq. (|B.ip into eg. (1^51) . 
eq. ([26| and eq. (|27|) can give the second-order rogue waves 
with five free parameters. Note that under the condi- 
tion J I » Jo, the second- rogue can split into three 
first-order rogue wave (triplets rogue wave) '6I'] rather 
than two. The dynamical evolution of ji^I^lp, jpl^lp 
and 77 for the parametric choice d — 1,6 = 2,ajo = 

-,_ft'o = 1, Jo = 0; Ji = 100 are plotted in the Figures 
7, 9 and 11 and their corresponding density plots are 
shown in the Figures 8, 10 and 12. There is another 
kind of second-order rogue wave, for example, jiJ^^lp is 
higher than second-rogue above. The dynamical evolu- 
tion of and 7y[^l for the parametric choice 

d = 2, & = 0, Wo = -, i^o = 1, >^o = 0, Ji = are plotted 

in the Figures 13-15. Note that eigenvalue Ai = A3 is 
real. The eigenvalue of rogue waves are different from 
the eigenvalue of solutions given in the past. 

Case 2). When n — 3, substituting eq. ljB.ip into ea. (l25t . 
eq. ([26l) and eg. ([ST)) can give the third-order rogue waves 
with six free parameters. Note that under the condition 
J2 >> Ji{i = 0, 1) or Ji >> Ji{i — 0, 2), the third-rogue 
can split into six first-order rogue wave rather. Circu- 
lar rogue wave [62| may be constructed by the condi- 
tion J2 >> Ji and J2 >> Jo. The dynamical evolu- 
tion of and fyt^' for the parametric choice 

d = 1,6 = 2,wo = -, JCo = 1, Jo = 0, Ji = 0, J2 = 8000 
are plotted in the Figures 16-18. At the same time, 
triplets rogue wave may be constructed by the condi- 
tion Ji >> J2 and Ji >> Jo. The dynamical evolu- 
tion of l^jl^lp and ryl^' for the parametric choice 

d = 1,6 = 2,a;o = ^,Ko - 1, Jo = 0, Ji ^ 100, J2 - 

are plotted in the Figures 19-21. Similarly, there is an- 
other kind third-order rogue wave, for example, 
is higher than third-rogue above. The dynamical evolu- 
tion of |£'['^lp, Ipl'^lp and ryl'^I for the parametric choice 

d = = 0,^0 = ^,i^o = 1, Jo = 0, Ji = 0, J2 = 
are plotted in the Figures 22-24. Note that eigenvalue 
Ai = A3 = A5 is real. The eigenvalue of rogue waves 
are different from the eigenvalue of solutions given in 
the past. According to analysis above, the n-order rogue 



waves may be controlled by n + 3 free parameters. 

From the above discussions, it is interesting to point 
out that the down amplitudes are dominant in the profile 
of rogue waves p^^l and f)[^l , so they are new kind of rogue 
waves when compared with the typical bright rogue wave 
E^^\ which are corresponding to the dark breathers in 
Figure 2 and Figure 3. So from our earlier understanding 
of breathers in other physical systems, we call these new 
type of solutions as dark rogue waves. Moreover, the dark 
rogue wave p'^' has one upper ring and three down peaks, 
and dark rogue wave fj^^^ has two lumps and one down 
peak. According to analysis, the n-order rogue waves 
must be generated by n-order breather solution. 

From the detailed literature on rogue waves, to the best 
of our knowledge, so far only bright rogue waves have 
been analyzed in detail but there is little report about 
the dark rogue waves in physical system. In the case of 
bright optical rogue waves, many results have actually 
connected the generation of supercontinuum generation 
(SCO) with rogue waves [63j. In recent years, the super- 
continuum white coherent source has attracted a lot of 
attention because of its potential applications in optical 
coherence tomography, spectroscopy, wavelength division 
multiplexing, etc. As reported in ref.[6J|, the modula- 
tional instability (MI) conditions for the generation of 
ultra-short pulses has already been investigated in the 
erbium doped nonlinear fibre and occurrences of noncon- 
ventional side bands have also been observed. This type 
of nonconventional side bands will be very useful to gen- 
erate large MI bandwidth which intern generates very 
short pulses. In this way, we believe that our rogue wave 
results in this paper can also be connected to the genera- 
tion of SCO. Similarly, the occurrence of dark rogue wave 
can also be connected to the results of ref. 64] , in the fol- 
lowing manner: In our above work 6J|, it has been shown 
that both bright and dark SIT solitons can be generated 
in the case of the anomalous and normal group velocity 
dispersion (GVD), in contrast to the well-established re- 
sults in the conventional fibre, where bright and dark soli- 
tons exists in the anomalous and normal GVD regions, 
respectively. From the above results, it is clear that the 
formation of dark rogue waves can also be connected in 
a similar way. Thus, it is interesting to analyze the rela- 
tion connecting the MI, SCG and rogue wave formation 
in optical system. 
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IV. CONCLUSION 

Thus, in this article, we have reported the rogue 
waves of the three physical fields E,p, and in a res- 
onant erbium-doped fibre system, which is governed by 
the NLS-MB equations. These rogue waves are con- 
structed by a Taylor series expansion of the correspond- 
ing breather solutions of the NLS-MB equations. As ex- 
pected, in contrast to the usual bright rogue wave E, we 
observe dark rogue waves for p and "q. The main feature of 
the dark rogue waves is the appearance of two (or more) 
dominant down peaks in its profile. In particular, there is 
one upper ring in the profile of the p, so it may be called 
as dark ring-rogue wave. The explicit form of 
and jyf"! are given by the determinant representation of 
the n-fold DT. The rogue waves in previous section can 
also be connected to the supercontinuum generation. 

As we have already described in the introduction, the 
singularity of the solutions generated by the DT is 
the main constraint to generate the dark rogue waves of 
the defocusing NLS equation. This perhaps shows that 
the dark rogue wave of the defocusing NLS equation can 
be investigated by other way such as by means of Hirota 
method. From the determinant representation of the T„, 
it is interesting to generate the higher order rogue waves 
so that the dynamical interactions of rogue waves can 
be analysed. 

In recent years, considering variable dispersion, vari- 
able nonlinearity and variable gain/loss, the investiga- 



tion of solitons in nonautonomous nonlinear evolution 
equation equations has also attracted a lot of attention 
among researchers [65l-l68l|. For example, Serkin and his 
coworkers have proposed a novel method to analyse the 
nonautonomous soliton equations [66,-£8j and the interac- 
tion of solitons in variable coefficient higher order NLS 
equation have been investigated in detailfos']. Using the 
results of the above papers and making use of our results 
in this paper, one can also construct the multi solitons, 
breathers and rogue wave solutions of the variable coef- 
ficient NLS-MB system. Moreover, it is also possible to 
obtain new type of rogue waves for other important cou- 
pled system in optics, such as the CH-MB equations [69| 
and variable coefficient CH-MB equations [70l Pzij . 
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V. APPENDICES 

Appendix I: In this appendix, we are giving 
expression for ti and elements of T„ 
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Appendix II: Determinant representation of Nth order DT is constructed in the form 
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with 
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Appendix III: Here, we are giving the construction of T of ^^.^(A) in detail. 
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Fig. 1: The first order breather ji^I^lp given with 

1 4 

specific parameters d — l,b — 2,u!o = -,ai = -. 

It o 

There are one upper peak and two caves in each 
periodic unit. 




Fig. 2: The first order dark breather IpI^'P for the 
values used in Figure 1. There are a upper ring and 
three down peaks in each periodic unit. 




Fig. 3: The first order dark breather ryl^l for the 
values used in Figure 1. There are two lumps and 
one down peak in each periodic unit. 




Fig. 4: The first order rogue wave jE'I-'^lp with 

specific parameters (i=l,6 = 2,ajo=2" There are 
one upper peak and two caves. 




Fig. 5: The first order dark rogue wave for 
the values used in Figure 4. There are one upper 
ring and three down peaks. 




Fig. 6: The first order dark rogue wave r/^^l for the 
values used in Figure 4. There are two lumps and 
one down peak. 



Fig. 7: The second order rogue wave given 
by eq. (P5|) with specific parameters 

d - 1, 6 - 2, wo - ^, i^o - 1, ^0 - 0, Ji - 100. 




Fig. 8: Contour plot of the wave amphtudes of 
I £'[2] 1 2 f-Q^. ^Yie values used in Figure 7. 




Fig. 11: The second order dark rogue wave Ty^^I 
given by ea. (|27p for the values used in Figure 7. 
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Fig. 12: Contour plot of the wave amplitudes of yyl^l 
for the values used in Figure 7. 
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Fig. 9: The second order dark rogue wave 
given by eg. ([25)) for the values used in Figure 7. 
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Fig. 10: Contour plot of the wave amplitudes of 
|p[2]|2 Jq^. ^Yie values used in Figure 7. 




Fig. 13: The second order rogue wave given 
by eq. (j25p with specific parameters 
1 



d = 2, 6 = 0, Wo = t:, ^0 = 1, ^0 = 0, Ji = 




Fig. 14: The second order dark rogue wave Ip'^'P 
given by ea.p6| for the values used in Figure 13. 



Fig. 15: The second order dark rogue wave 77^' 
given by ea. (l?fl) for the values used in Figure 13. 




Fig. 16: The third order rogue wave jijl^lp given by 
eq. (|25p with specific parameters d — l,b — 2,ll!q = 

l,Ko^l,Ja = 0,Ji= 0, J2 8000. 



Fig. 19: The third order rogue wave jE'I'^lp given by 
eg. (PS)) with specific parameters d = 1, 6 = 2, = 

^,Ko = 1, Jo = 0, Ji = 100, J2 = 0. 




Fig. 20: The third order dark rogue wave 
given by eg. (1261) for the values used in Figure 19. 




Fig. 17: The third order dark rogue wave jpl'^lp 
given by eg.([26| for the values used in Figure 16. 




Fig. 18: The third order dark rogue wave Ty^'^I given 
by eq.(|27)l for the values used in Figure 16. 




Fig. 21: The third order dark rogue wave ry^'^I given 
by eg. ([27]) for the values used in Figure 19. 




Fig. 22: The third order rogue wave given by 

eq.([25|)with specific parameters 

d=^,b^0,ujn^l,Ko = l,Jo^0,Ji^0,J2 = 0. 




Fig. 24: The third order dark rogue wave yy^'^I given 
by eq.(|27p for the values used in Figure 22. 



